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Abstract: Traditional methods are difficult for the measurement of refractory metals in high temperatures. The 


containerless electrostatic levitation method is combined with image measurement technique. During the heating 


process, the levitated sample images are recorded by CCD camera. The image edge is extracted by direction 


gradient algorithm, and the volume can be calculated by a spherical harmonic function. Then thermal expansion 


coefficients at different temperatures can be calculated. The pure Zirconium sample is used to test the measurement 


of thermal expansion coefficient in the range of 750-1800°C, the measuring method and accuracy are verified. 
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0 引言 样品 ， 再 利用 光学 等 方法 测量 样品 在 不 同 温 度 下 
胀 系数 是 材料 基本 的 热 物 性 参数 之 一 ， 是 。 ”长 度 或 体积 的 变化 。 由 于 测量 过 程 样品 与 测试 仪 
工程 设计 、 精 密 仪器 制造 、 焊 接 、 材 料 制 和 等 领 。 ” 涡 的 接触 ， 能 加 热 到 的 最 融 温 度 有 限 ， 并 且 在 温 
域 中 的 不 可 忽略 的 重要 参数 了。 度 较 高 时 材料 样品 的 温度 分 布 均匀 性 差 ， 影响 测 
传统 的 热膨胀 系数 的 测量 方法 包括 项 杆 法 、 干 。 ”时 精度 ， 导 致 有 效 的 温度 测量 范围 窄 (一 般 温度 
涉 法 和 热膨胀 法 等 中 ， 这 些 方法 都 是 采用 加 热 装 测量 范围 低 于 1500K) 。20 世 纪 60 年 代 ， 以 电磁 县 
置 与 测量 材料 样品 接触 通过 热传导 的 方式 加 热 浮 技术 为 代表 的 无 容器 非 接触 测量 技术 开始 用 于 
材料 热 物 性 的 测量 。 无 容器 技术 利用 气流 ”、 声 
基金 项 目 ， 国 家 国际 科技 合作 专项 2011DFA10440)， 中 国 场 ”、 电 磁场 ”、 静 电场 ”等 手段 使 样品 悬浮 起 
200028 来 ， 脱 离 与 容器 壁 的 接触 ， 通 过 大 功率 的 激光 器 
生 ， 专 业 方向 为 静电 悬浮 无 容器 材料 实验 技术 。 熔点 使 样品 熔化 ， 另 外 巧 浮 的 样品 不 易 受 器 壁 的 
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迹 挡 ， 通 过 图 像 技术 能 够 更 高 精度 地 测量 样品 体 
积 的 变化 。 对 于 温度 范围 绞 高 的 金属 材料 的 热 膨 
胀 系数 测量 ， 需 要 在 高 真空 环境 下 ， 一 是 真空 环 
境 下 样品 只 通过 辐射 散热 可 以 加 热 到 更 高 的 温 
度 ， 二 是 防止 样品 氧化 影响 测量 精度 。 气 悬浮 和 
声 臣 浮 不 适用 于 真空 的 环境 ， 电 磁 惹 浮 技 术 的 其 
浮 和 加 热 是 相 耦 合 的 ， 不 能 独立 控制 加 热 的 温度 ， 
对 于 低温 范围 无 法 测量 ， 并 且 电 磁场 产生 的 涡流 
式 使 熔化 的 样品 形状 不 稳定 ， 影 响 测量 精度 。 
本 文 利用 静电 基 浮 无 容器 材料 实验 装置 "进行 
热膨胀 系数 的 测量 。 静 电 甚 浮 技术 解决 电磁 悬浮 
带 来 的 问题 ， 静 电 巧 浮 的 加 热 是 独立 控制 的 ， 样 
品 从 加 热 到 熔化 都 处 于 稳定 的 状态 ， 该 装置 使 用 
CCD 相 机 与 背景 光源 , 从 侧面 获取 清晰 的 样品 图 像 
并 通过 图 像 处 理 算法 计算 样品 的 体积 变化 。 最 后 ， 
通过 对 纯 错 样品 的 实验 ， 对 该 测量 系统 进行 了 验 
证 和 测试 。 

实验 装置 

对 材料 热膨胀 系数 的 测量 是 基于 静电 悬浮 材料 
科学 实验 平台 ， 该 实验 装置 包括 真空 腔 体 、 悬 浮 
位 置 控制 系统 、 加 热 激光 器 、 双 波长 测 温 仪 和 CCD 
测 温 仪 ， 装 置 的 结构 如 图 1 所 示 。 真 空 腔 体 通过 一 
组 机 械 泵 和 分 子 泵 保持 在 5X10 Pa 以 下 的 高 真空 
的 实验 环境 。 测 试 的 材料 为 直径 2-3mm 的 球形 样 
品 ， 在 真空 腔 体 内 部 有 一 组 间距 为 10mm 的 高 压 平 
行 极 板 ， 产 生 的 高 压 电 场 使 带电 荷 的 测量 样品 在 


库仑 力 的 作用 下 , 抵消 样品 自身 的 重力 悬 译 起 来 ， 
样品 的 位 置 通过 反馈 控制 系统 保持 稳定 ”。 


一 一 


CCD 相机 nie 


图 1 静电 悬浮 装置 结构 图 


样品 稳定 悬浮 之 后 ， 通 过 激光 加 热 器 对 样品 进 
行 加 热 和 双 波 长 测 温 仪 获取 样品 的 温度 。 对 测量 


样品 的 加 热 采 用 808nm 的 大 功率 半导体 激光 器 , 激 
光 器 的 功率 为 90W， 激光 通 过 光纤 引出 ,经 过 透镜 
聚焦 到 样品 表面 ， 对 于 金属 材料 最 高 加 热 温 度 可 
以 达到 2500C 。 样 品 加 热 温度 的 测量 通过 一 台 采 
样 频率 为 50Hz 的 双 波 长 测 温 仪 〈SENSORTHERM 
MQ22) 实现 。 双 波长 测 温 仪 的 利用 两 个 临近 波长 
的 辐射 度 的 比值 计算 样品 的 温度 ， 所 以 不 与 样品 
的 发 射 率 、 测 量 立体 角 等 参数 有 关 ， 测 温 精 度 较 
高 。 材 料 的 热膨胀 系数 是 在 加 热 的 过 程 中 通过 图 
像 采集 装置 进行 测量 。 


图 2 CCD 相 机 采集 到 的 加 热 样品 的 图 像 


当 样 品 在 加 热 到 很 高 温度 之 后 ， 自 身 就 会 产生 
明亮 的 可 见 光 ， 使 图 像 边 缘 变 得 模糊 。 为 了 获得 
清晰 的 图 像 轮廓 ， 采 用 上 臣 素 灯 和 毛 玻璃 作为 背景 
光源 ， 减 弱 样品 自身 发 光 的 影响 。 由 于 样品 处 于 
无 容器 悬 译 状态 ， 即 使 样品 熔化 为 液态 ， 在 表面 
张力 的 作用 下 仍 保持 近似 球形 ， 根 据 球形 轴 对 称 
的 性 质 ， 采 集 侧 视 的 图 像 就 可 以 通过 图 像 处 理 算 
法 精确 得 到 样品 的 体积 。 


2 图 像 处 理 算法 


为 了 根据 样品 的 图 像 精确 测量 样品 的 体积 ， 需 
要 通过 图 像 处 理 的 方法 实现 。 首 先 ， 要 从 背景 中 
提取 样品 的 边缘 ， 边 缘 提 取 根 据 像素 点 灰 度 值 的 
梯度 变化 来 确定 边缘 位 置 “， 但 梯度 是 具有 方向 
的 ， 通 用 的 方法 采用 水 平和 垂直 两 个 正 交 的 方向 
计算 梯度 变化 。 由 于 样品 为 球形 ， 在 样品 图 像 中 
垂直 于 边缘 各 点 切线 的 方向 具有 最 明显 的 梯度 变 
化 ， 所 以 只 采用 水 平 或 垂直 方向 进行 梯度 变化 的 
检测 ， 对 边缘 的 提取 精度 不 够 ， 本 系统 针对 圆 的 
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边缘 特点 采用 了 方向 梯度 算法 进行 边缘 像素 点 提 
取 。 

采用 360 组 从 圆心 Co, 加 ) 发 出 的 方向 向 量 , 每 一 
组 向 量 方向 上 将 采集 到 一 组 具有 方向 特性 的 像素 
序列 。(%, yo) 是 用 来 产生 方向 向 量 的 近似 圆心 位 
置 ， 可 以 采用 外 切 和 矩形 框 的 中 心 近似 蔡 代 。 产 生 
的 360 组 方向 向 量 上 的 像素 序列 表示 如 下 : 


Xi =7,sinQ + Xo 


yj =7;COsO + yor,i=1,...,360 (1) 
0.=2xi/360 


式 中 i 是 360 组 方向 向 量 的 索引 ，j 是 每 一 个 方 
向 向 量 上 像素 序列 的 索引 , x 和 7 是 水 平和 垂直 的 
像素 点 坐标 , 7 是 产生 像素 序列 的 半径 , 由 于 图 像 
的 像素 点 是 离散 的 ， 半 径 六 以 一 个 像素 点 进行 递 
增 ， 直 到 产生 的 像素 点 到 达 图 像 的 边缘 。 

之 后 分 别 对 每 一 组 两 维 像素 点 序列 作为 一 维 的 
像素 灰 度 值 进行 处 理 ， 对 像素 的 灰 度 值 进行 微分 
得 到 该 方向 上 梯度 变化 ， 样 品 的 边缘 就 是 在 每 一 
组 方向 向 量 上 最 大 灰 度 梯度 变化 的 像素 点 位 置 ， 
处 理 得 到 360 个 边缘 像素 点 坐标 (5,y) 。 

为 了 提高 边缘 提取 的 精度 , x 和 ?方向 坐标 的 平 
均值 用 来 作为 新 的 近似 中 心 点 (6,yo)， 用 这 个 新 
的 中 心 点 重复 式 〈1) 的 过 程 ， 重 新 产生 360 组 向 
量 并 搜索 边缘 像素 点 。 最 终 应 用 方向 梯度 算法 得 
到 样品 图 像 边缘 提取 的 结果 如 图 3 所 示 。 
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图 3 方向 梯度 算法 提取 样品 图 像 的 边缘 


通过 方向 梯度 方法 得 到 边缘 信息 为 一 组 像素 点 
坐标 ， 这 些 像 素 点 坐标 是 离散 的 ， 即 使 采用 高 分 
状 率 的 CCD 相 机 进行 采集 , 像素 点 离散 的 特性 也 无 
法 完全 消除 提取 的 图 像 边缘 的 锯齿 状 ， 这 对 精确 
计算 样品 的 体积 产生 影响 。 

系统 采用 球 谐 函数 拟 合 的 方法 ， 通 过 有 限 的 离 
散 像素 坐标 点 ， 通 过 拟 合 得 到 精确 的 球 谐 函 数 表 
示 。 球 谐 函 数 的 表示 形式 如 下 : 


6 
R(0) = 2c,P" (cosO) (2) 

n=0 
式 中 ，Pr"(cos9) 是 阶 数 为 n 的 Legendre 多 项 式 ， 
阶 数 越 高 拟 合 的 精确 度 越 高 ， 但 计算 量 也 随 之 增 
大 ， 系 统 采 用 6 阶 进行 拟 合 ，c, 是 需要 拟 合 过 程 需 

要 计算 的 系数 。 

为 了 对 球 谐 函数 进行 拟 合 ， 需 要 把 边缘 提取 得 
到 的 像素 点 从 直角 坐标 转换 为 如 下 的 极 坐标 形 


式 : 
R= -xX) P+, -hPL 
站 


i 


式 中 ， 到 和 到 是 水 平和 垂直 方向 上 的 像素 度量 
系数 〈 训 米 /像素 ) ， 表 示 图 像 两 个 像素 点 之 间 的 
实际 距离 ， 和 相机 CCD 蕊 片 、 与 样品 的 距离 以 及 镜 
头 焦距 等 参数 有 关 ， 通 过 对 标准 尺寸 的 样品 进行 
标定 得 到 ，Xo 和 部 是 精确 的 像素 坐标 中 心 ， 需 要 
通过 拟 合 确定 。 由 式 (3) 可 以 把 样品 图 像 边缘 点 
由 直径 坐标 (%,y) 转化 为 极 坐标 (及 ,2) 来 表示 。 

通过 式 (3) 表示 的 极 坐标 点 对 式 〈2) 中 的 系 
数 进行 拟 合 ， 拟 合 系数 包括 Legendre 多 项 式 的 6 
阶 系数 co 到 cs 以 及 坐标 中 心 X 和 成 » 可 以 通过 对 
如 下 的 目标 函数 求 取 最 小 值 来 得 到 相关 的 拟 合 系 
数 : 


360 2 


(XU 马 ,coco)=>[R-RO) 0) 
iL 


式 中 ， 在 角度 2 的 方向 上 ， 样 品 图 像 边缘 点 与 
对 应 的 球 谐 函数 的 半径 差 值 进行 平方 , 对 360 组 数 
据 进 行 累加 ， 对 应 目标 函数 的 最 小 值 就 是 拟 合 系 
数 的 最 优 解 。 首 先 ， 坐 标 中 心 X 和 郊 用 各 像素 点 
水 平和 垂直 坐标 的 平均 值 进行 初始 化 ， 根 据 式 
(2) ，R(O) 对 于 拟 合 系数 是 到 cs 是 线性 的 ， 可 
以 采用 最 小 二 乘法 进行 计算 ， 得 到 初始 Xo 和 郊 下 
的 最 优 系数 解 。 分 别 改变 Xo 和 页 的 值 ， 每 次 改变 
都 需要 重新 对 式 (3) 中 的 极 坐标 点 进行 计算 ， 并 
采用 最 小 二 乘法 重新 计算 使 目标 函数 最 小 的 一 组 
新 的 Legendre 多 项 式 系数 。 经 过 拟 合 得 到 球 谐 函 
数 式 (2) 的 最 优 系数 之 后 ， 可 通过 如 下 公式 对 样 
品 图 像 轮廓 半径 的 积分 计算 得 到 样品 的 体积 : 


过 | 
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加 热 激光 功率 (W) 
口 虽 寺村 名品 


i | 
V 本 ， R(O) sinGd0 (5) 
得 到 样品 的 体积 之 后 ， 结 合 双 波 长 测 温 仪 采集 
温度 变化 数据 ， 可 以 通过 如 下 公式 计算 样品 材料 
的 膨胀 系数 : 
1 ,AV 


pb = a (6) 


3 ”实验 过 程 和 结 


在 实验 开始 之 前 ， 首先 悬浮 一 个 标准 尺寸 的 不 
锈 钢 球 样品 对 图 像 采集 系统 的 分 辨 率 进行 标定 ， 
标定 的 样品 莽 浮 稳定 之 后 , 调整 CCD 相 机 镜头 的 光 
图 、 焦 距 和 图 像 的 对 比 度 等 参数 ， 使 样品 图 像 清 
晰 ， 一 旦 调整 好 在 整个 实验 过 程 中 都 保持 这 些 光 
学 参数 固定 不 变 。 通 过 测量 标准 样品 在 图 像 中 的 
像素 数 ， 得 到 式 (3) 中 水 平和 垂直 方向 上 的 像素 
度量 系数 和 。 

系统 经 过 标定 后 ， 放 入 经 过 真空 电弧 炉 制备 的 
纯 错 样 品 ( 纯 度 99. 99%) ， 样 品 稳定 悬浮 状态 下 ， 
通过 半导体 激光 器 加 热 直到 完全 熔化 ， 然 后 关闭 
加 热 激 光 ， 样品 通 过 辐射 降温 ， 紧 接着 凝固 再 辉 。 
测 温 仪 得 到 的 熔点 温度 (T=1852'C ) 与 钳 的 熔点 
相符 合 。 实 验 过 程 的 激光 加 热 和 测 温 仪 测量 得 到 
的 温度 数据 如 图 4 所 示 。 
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图 4 实验 过 程 样品 加 热 和 测量 温度 数据 曲线 


测 温 仪 温度 的 采集 周期 为 20ms， 系 统 每 次 从 测 
温 仪 获取 一 次 温度 数据 ， 就 同步 触发 一 次 图 像 采 
集 ， 这 些 图 像 数 据 连 同 对 应 的 温度 信息 一 起 存储 
在 计算 机 的 硬盘 上 ， 待 实验 结束 后 对 图 像 进行 处 
理 。 由 于 装置 的 悬浮 样品 的 位 置 稳定 精度 在 士 
0. 2mm, 可 以 根据 第 一 幅 图 像 选 定 样品 图 像 的 处 理 


体积 (mm’) 


区 域 并 通过 批 处 理 的 方法 对 每 幅 图 像 进 行 处 理 ， 
根据 式 〈5) 计算 样品 的 体积 ， 结 合 温度 信息 得 到 
样品 体积 随 温 度 变 化 的 曲线 ， 如 图 5 所 示 。 
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图 5 样品 体积 随 温度 变化 曲线 


在 样品 体积 测量 过 程 中 ， 低 温 阶段 和 熔点 温度 
附近 测量 精度 较 差 ， 不 对 其 进行 计算 ， 材 料 的 热 
膨胀 系数 根据 式 6) 和 体积 与 温度 的 测量 数据 进 
行 计算 ， 得 到 样品 的 热膨胀 系数 与 温度 的 关系 如 
图 6 所 示 。 
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图 6 样品 热膨胀 系数 的 变化 曲线 


由 于 热膨胀 系数 的 计算 与 样品 的 质量 参数 无 
关 ， 所 以 其 精度 只 取决 于 图 像 测量 样品 体积 和 双 
波长 测 温 仪 的 精度 。 利用 最 小 二 乘法 对 图 6 中 的 数 
据 进 行 线性 拟 合 ， 得 到 固态 样品 铬 的 热膨胀 系数 
随 温度 变化 的 表达 式 为 : 

Ps, =(1.573x10°)—(3.843x10”)xT+2.3% 

(750<T<1800°C) 

式 中 ， 样 品 热 膨胀 系数 局. 的 单位 为 K"， 样 品 

温度 7 的 单位 为 'C 。 
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4 结束 语 
本 文通 过 静电 悬浮 无 容器 的 方法 测量 金属 材料 
的 热膨胀 系数 ， 提 高 了 传统 测量 的 方法 的 温度 测 
量 范 围 ,最 高 有 效 温度 达到 2000K 以 上 ,并 且 比 电 
磁悬浮 等 其 他 无 容器 测量 方法 精度 更 高 ， 通 过 背 
景 光 源 提高 采集 样品 图 像 边缘 的 对 比 度 ， 利 用 基 
于 梯度 的 边缘 提取 算法 和 球 谐 函 数 拟 合算 法 进 一 
步 提 高 样品 体积 的 测量 精度 。 最 后 ， 通 过 金属 钳 
样品 的 实验 验证 了 测量 方法 的 有 效 性 ， 为 各 领域 
的 研究 人 员 提 供 精 确 的 高 熔点 合金 材料 的 热膨胀 
系数 测量 数据 和 方法 。 
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